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LITERAL EXPRESSION FOR THE MOTION OF THE MOON'S 

PERIGEE. 

By De. G. W. Hill, West Nyack, N. Y. 

The earlier investigators of the Lunar Theory contented themselves with 
giving numerical values for this quantity. The results of Clairaut, D'Alembert, 
Euler, Laplace and Damoiseau are of this nature. Beyond the rudest approxi- 
mation, Plana was the first to give the value in a literal form. This was nearly 
reproduced by Pont^coulant. The salient portion, which is a function of the 
ratio of the month to the year, is given by these two authors to the order of 
the seventh power of this ratio. Delaunay, in his treatment of the subject, 
(Oomptes Rendus, t. LXXIV, p. 19) has added the two following terms, 
viz. those which involve the eighth and ninth powers of the mentioned ratio. 
The correctness of these has, however, been called in question by M. Andoyer 
(Annates de la Faculte des Sciences de Toulouse, t. VI) who has also given 
other values for them. They may be seen in Tisserand's Mecanique Celeste, 
t. Ill, p. 412. I have not been able to consult M. Andoyer's memoir, and 
do not know what method he used in obtaining his results. The comparison 
I made at the end of my memoir " On the part of the motion of the Lunar 
Perigee, etc. (Acta Math., Vol. VIII) of Delaunay 's series with my numerical 
value, indicated with some probability, that the newly added terms were, one 
or both, too large ; which corresponds with what M. Andoyer has found. In 
this state of matters I have thought it would not be without interest to test the 
validity of M. Andoyer's corrections ; and I have determined to add two more 
terms to the series, viz. those factored by m 10 and m 11 . 

For this purpose, I shall employ the method of my above-mentioned 
memoir. There the computation was carried out in the numerical fashion, 
here it is proposed to give algebraic developments. It is there shown that the 
determination of the lunar inequalities of the type of the evection and the 
motion of the perigee depend on, at least as far as a first approximation is 
concerned, the integration of the linear differential equation of the second 
order 

D 2 w = <9w , 

where w is the unknown, 6 a periodic function of double the mean angular 
distance of the Moon from the Sun, involving only cosines, and D is an 
operator such that D (a£') = v«C, C being the trigonometrical exponential 
corresponding to the mentioned mean angular distance. The motion of the 
perigee depends solely on the coefficients of 8, and these can be found when 
Annals of Mathematics. Vol. IX, No. 2. 
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we know the coefficients of the inequalities of the type of the variation. I 
have given the latter in a literal form to the 9th order inclusive (Anier. Jour. 
Math., Vol. I, pp. 142-143). 

We adopt a moving system of rectangular coordinates, the origin being 
at the centre of the Earth, the axis of x constantly passing through the centre 
of the Sun, and, in place of x and y, we use the imaginary coordinates 

M=a; + yi/' — 1, 8 = x — yy y —\. 
Then 

u=£a t ^+ l , s=2'a i C" 2<_1 . 

where, in the summation, i receives all integral values from — oo to -\- oo , zero 
included, and the a 4 are constants, being each equivalent to the same constant 
multiplied by a function of the ratio of the month to the year, which is of the 
|2i|th order with respect to this parameter. For simplicity in writing, then, 
we assume that the value of a is unity ; consequently, as written, a t always 

denotes — . 

a 

In pushing the development of 6 to the degree of approximation we de- 
sire, the values of the & t given (Amer. Jour. Math., Vol. I, pp. 142-143) gen- 
erally suffice ; but it will be perceived from the approximate expression for 8 
(Motion of Lunar Perigee, p. 13) that it will be necessary for the determination 
of c, the ratio of the anomalistic to the synodic month, to the 11th order 
inclusive, that we should know the term factored by m 10 in a, +-a_j ; it is not 
necessary, however, that ai and a_! separately should be known to this degree 
of approximation. Hence, we now proceed to obtain this term. From the 
equations given (Am. Jour. Math., Vol. I, p. 137) by neglecting all terms whose 
order exceeds the 10th we derive 

3 (2 + m) m 2 ,, . n N 3(1 — m) m 2 , , . „ . n v 

a, -f a_, = — - v , — '- — , (1 + 2 a,a_i) — . v , '- — , (a_i + 2a_ 2 + 2a,a_ 3 ) 

11 l 6 — 4m + m 2 v ' ' 6 — 4m + m 2 v 1 2 l 3 ' 

22 — 4m + m 2 , . v Q 

~ 6 - 4m + m 2 ^^ + ^^ ~ 9 aA ' 

By substituting in the right member of this the values of the & t in powers of 
m it is found that the term in m 10 in a, + a -i is 



, 1605921808447 10 
+ 2 18 . 3 8 . 5 3 m 
We have 



« = -\jz + m 2 ] + 2[i) log ^p g + m]- [D log v'DuBsf 



D 2 [log v'DuDs] 
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The first term of this can be developed from the formula 

x , , IPu -\- 2mDu -4- 4 m 2 s . . , i , o r s 2 . v e> *2i 
-3 + m 2 = !— » + |m 2 = l + 2m + 4jm 2 +2. R,^' , 

I It 

where we have R_ t — R t . The equations which determine the values of the 
R t , i? with an error of the 12th order, Ii l with one of the 11th order, R 2 with 
one of the 9th order, and R 3 with one of the 7th order, are 

#0 + ( a i + a -i) -#1 + ( a 2 + a_ 2 ) R 2 = I m 2 a_, , 
a_!i? + (1 + a_ 2 ) R x + (a! + a_ 3 ) R 2 + a 2 i? 3 = — 4ma_, + fm 2 , 

a_ 2 -ff + (a_ x + a_ 3 ) R t + R z + a. x R s = 8 (1 — m) a_ 2 + f mX , 

a-1^2 +-Rs= 24 »-s + f m 2 a 2 . 

Solving these by successive approximations in using the known literal values 
of the a 4 , we get 



%—•-■+ 4* + ¥-■+ 1^+ 2 -^»- + j^-f 



, 1260881 10 291394307 u 

J- m " — • m 

^ 2 9 .3\5 2 7 .3 6 .5 3 ' 

r> 3 2 . 19 3 . 20 4 , 43 , , 18709 , , 759413 7 
-^ 1 = g m + W + T + W + 2T3 3 " + 2"».3 4 .5 

, 6675059 8 41991161 , 4528083484913 w 



2 8 .3 5 .5 2 2 7 .3 6 .5 3 2 17 .3 7 .5 4 

2937 5 , 23051 s , 97051 , , 

2 6 .5 ^2 4 .3.5 2 ^ 2 s . 5 3 ^ 2 10 .3 3 .5 4 



„ 33 4 , 2937 5 , 23051 - , 97051 - , 334413271 8 



73 xdyo * 
^3 = -gj- m 6 • 

We next attend to the coefficients of -jr~ — Z . U i £ 2< . The formulas given for 

these (Motion of Lunar Perigee, p. 12) in general suffice ; it is necessary, how- 
ever, to push the development of U x and U_ x farther, so that terms of the 10th 
order may be included. Thus their equivalents read 

U x = 2 [h, — h_, h 2 + V h_, + 2h, 3 h_, 2 — h, 3 h_ 2 — 3h, h_, 2 h 2 + 2h, h 2 h_ 2 

+ h_ 1 2 h 3 -h_ 2 h 3 ], 

U_ x = — 2 [h_, — h, h_ 2 + h_! 2 h, + 2h_, 3 h L 2 — h.! 3 h 2 — 3h_! h t 2 h_ 2 

+ 2h_ 1 h_ 2 h 2 + h 1 2 h_ 3 -h 2 h_ 3 ]. 
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By substituting in these equations the values of the a ( in powers of m, and 
making the assumption that the coefficient of m 10 in a, is 0, which cannot lead 
us into error within the limits we set to the approximation, and putting 

we get 

. 5 2 1 s . 5 4 , 43 5 318575 „ 2297593 7 

A = - a*™ 2 - 2T3 m + r* m + F73* m - 1F73T m ~ <F7F75 m 

9225887 8 3471983789 9 12903700736069 10 
m 8 na n , p. m 9 — „,„ „,, „ A m 10 , 



2 5 . 3 6 . 5 2 2 9 . 3 7 . 5 3 2 19 . 3 8 . 5 4 

n _ 7 . 19 3 , 53 < 155 . 12941 , 904921 . 
Bx — 2~ 2 m + 273 m + 273~ 2 m + FTP m _ W7W m _ 2 9 TF75 m 

35308207 8 2190838913 9 . 29589760583167 I0 
2*. 3 6 . 5 2 2 10 . 3 7 . 5 3 ^ 2 18 . 3 8 . 5 4 ' 

>! - 265 m < 4. 1067 m * 4 38261 m a , 755591 , 405840581 . 
2 — 2 7 2 5 5 2 4 3 2 5 2 2 6 3 2 5 3 "2 15 3* 5^ ' 

„ 3 4 403 5 3773 6 246139 7 1077852389 8 

^2— 2 2™ 2 4 .3.5 2 3 .3 2 .5 2 2 5 .3 3 .5 3 2 12 .3 4 .5 4 ' 

A - 1677 m* 

r> 2431 6 

^s = -210- m • 

The coefficients of the function 8 have the following equivalents : 

6 = 1 + 2m ~\ m 2 + 4(A 2 + -4/) + 2 (^ + B*)-R* , 
8, = 4 (1 + m) A, - 2B t + 4 (^^ 2 + 4,4,) + 2 (^ + ^5.) - R, , 
2 = 4 (1 + m) 4 2 - ±B 2 + 2 (4 X 2 + 24^) - B* + 2^, - R 2 , 
6> 3 = 44 3 - 6B S + 44 t 4 2 — 2B X B 2 - R 3 . 

When the expressions in powers of m for the quantities A, B, and R are sub- 
stituted in the preceding equations, the results are 

a i , o i 2 , 255 4 . 1Q 5 . 80 6 . 533 7 , 11230225 8 

#o = 1 + 2 m — i m + -gr m 4 + 19m 5 + -g- m 6 + g-g, m 7 + ^ g 3 m 8 

, 1576037 „ , 49359583 , . 720508007 u 
+ -2 7 T3^ m +-2T3 5 - m +-2TF7F m ' 
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a 15 , 57 , „ ( 23 , 68803 „ 1792417 7 
01 = ~ T m "2^ m 11m 4 — 273 m — 297F m 2 10 . 3 s 

7172183 8 596404499 , 2641291011773 10 
2 7 . 3 4 . 5 2 9 . 3 s . 5 2 2" . 3 6 . 5 3 ' 

a 111 4 , 1397 , . 8807 6 , 319003 7 . 252382507 8 
«.= -gr ™ 4 + -gr- m * + 2^7375 m 6 + g^-^-g, m? + 2 ,„ yy m», 

a H669 , 

We employ now the system of equations given (Motion of the Lunar 
Perigee, p. 14) and for brevity of notation put [i] for (c -\- if — 6 . The 
equations, written to the requisite degree of approximation, are 



[- 3] b_ 3 - # x b_ 2 - 6 2 b 

- 0,b_, + [- 2] b_ 2 - # t b_ 

- 2 b_ 3 - 6,b_ 2 + [- 1] b_ 

- 6» s b_ 3 - 2 b_ 2 - ^b. 

- 3 b_ 2 - 6 2 b 



- 8 3 \ = o : 



- # 2 b - e B b 1 = o , 

- 0,b, - # 2 b x - e 3 b 2 = o : 

+ [0] b - fl.b, - # 2 b 2 = : 

- ^b + [1] b t - ^b 2 = : 

- 6 2 b - Bfr + [2] b 2 = , 



That the relative degree of importance of the terms of these equations may be 
perceived, it may be pointed out that the diagonal line of coefficients [ — 3], 
[ — 2], . . . , [1], [2] are all of the zero order of magnitude except [ — 1] and 
[0], the first of which is of the first order, and the second of the third order ; 
but the latter we need not concern ourselves about ; and S t is of the 2ith 
order. From this it follows that, if we write the quantities b in the order b_ u 
b[, b_ 2 , b 2) b_ s , leaving out b , which is an arbitrary quantity, the first is of the 
first order, and every succeeding one an order higher, so that b._ 3 is of the fifth 
order. In order to have the equation determining c, it is necessary to elimi- 
nate the 5 mentioned b's from the group of equations. The readiest method 
of accomplishing this is to proceed by successive approximations using for- 
mulas of recursion. To attain the desired degree of accuracy, three approxi- 
mations are necessary, each of which will give three terms in powers of m in 
the value of each b involved. When the values of the five b's have been obtained 
and substituted in the middle equation, after the rejection of the useless factor 
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b , we have the following equation serving for the determination of c : 

— %[ u ± i_ t 1 0*0 

L[- 1] [1] t- 1] [- 2] + [1] [2] J Wl * 2 

~ LF 7 ^ + PI J 6i ~ [=m[- 1] C- 2] 2 + [- 2] 2 [- 3]] ^ 



_,Q 



+ r ,-,, r ;-,r on + 



[- l] 2 [- 2] 2 ' [- l] 2 [1] [- 2] ' [- 1] [l] 2 [2] 



+ r ii2 r on r Qi + r n r _ gp [_ 3] J ^^ 



1 

[_!]»[_ 2] [-3] ' [-1][- 



r l . i_ l . 

I r \~\ [ 0.12 "+' r 112 rn "' r _ i i riu "i 



[-2]* ' [-1] 2 [1] ' [-1J[1] 2 ' [-l] 2 [-3] 



2 2 2 

[- 1] [1] [- 2] + [- 1] [1] [2] + [- 1] [- 2] I 



~ t-^tt-^ + r^™ + r-^n^Vpnyj] W 



2 L[- 1] [- 2] [- 3] + [- 1] [1] [- 2] + [- 1] [1] [2]] W 

- 2 [f^ + FTTF^]^^°- 

From this all terms unnecessary to the desired degree of approximation have 
been rejected. 

It appears desirable to give some details as to the treatment of the fore- 
going equation. First we form the various products of the 8 involved ; each 
is limited to the terms needed for the degree of approximation wished. 



225 , ,855 . , 5889 .,.,„, , , 697679 „ , 853817 



0, _ -^ m +-22-m + - w ~ m + d71m + -gr-g- m + -gs-g 

235899233 10 1733519201 „ 19979134939549 12 
+ 2 n ; 3 s m + 2 9 . 3 4 . 5 m H 2 18 . 3 5 . 5 2 m ' 

ai 50625 8 . 192375 9 . 2787075 10 , 6370695 „ 

e i = — §r- m + 2 s + 2 s — + — 2 s — 

. 353456169 12 . 649258747 18 

"T ffio m H SI5 m . 



m 
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at 11390625 u , 129853125 13 , 2868159375 M 
w i = 2 s m + 2 s m 2 s m ' 

ft2ft __ 24975 „ 693945 , 1188267 10 21446525 „ 

°i u i g 8- ~*~ — 2 s — 2* — 2 15 

, 4710472379 12 
+ 2».8.5 m ' 

„, _ 12321 . 155067 , 2018553 3 10 85123117 u 
2 2^ H 2 1- + 2 12 .5~ + 2 r r875 r ' 

„,„ _ 6619375 u 241552125 13 g lffl _ 2772225 M 55958985 „„ 

"i ^2 gs — "•" 2 15 ' 2 2 15 — 2" ' 

e » g _ 39382875 12 6 60 - 19428885 m « 

In the next place, by neglecting quantities of the 7th order, the equation 
may be written 

and, if we put 

it can be given the form 

[O] 2 + 2(0 O -1)[O] + #? = 0, 
whence is derived 



c 2 = l + i/(^ o -l) 2 -0' 1 2 . 
By substituting the previously given developments of O and X we get 

, , , o 1 2 225 3 3135 t 139973 , 4611319 , 
c 2 = 1 + 2m - g m "2T- m %r 2 15- FTS" ' 

. 3 2 201 3 2367 4 111749 „ 4095991 „ 

o = 1 + m - ^ m 2 - -g- m 3 - -^ m* - -^ m* - -^-3- m« . 

These equations are correct to the last power of m set down. 

This value of c may be substituted in all the terms but the two first of the 
equation which determines it ; and the latter is thereby reduced to a manage- 
able form. The values of the reciprocals of the quantities denoted by the 
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symbols [ — 1], [1], [ — 2], [2], [ — 3], developed to the needed degree of 
approximation, are 

1 1 , 3 213 2259 2 70973 s 3501259 . 
= — n-rn- 1 — 7U — -^j-m — - 7im - m 2 =^— m 3 „,« n m 4 , 



[_ 1] — 2 2 2 4 2 8 2 10 2 13 2 15 . 3 

1 1 1 , 5 2 , 563 s , 6119 4 

1 1,1 , 1 2 643 3 10807 t 532047 , 
- T3 + S3 m + as m — w m an- m sis— m . 



[— 2] ~~ 2 3 ' 2 3 ' 2* 2 10 2 

i i i , 13 2 . 8557 

™ + S5-55 m + HH-oi m ' 



[2] 2 3 . 3 2 s . 3 2 ' 2 5 . 3 3 ' 2 10 . 3 4 
1 1,1 

Q3 Q ~ 04 Q m 



[—3] 2 3 .3 ' 2 4 .3 

The substitution being made in the eight terms of the left-hand member of the 
equation, the result follows, in which, for facility of verification, we write each 
fraction separately and in the order in which it arises from each of the eight 
terms. 

50625 6 , r 1022625 , 24975 



F ml + 



r 1022625 , 24975T T 

L — t*— + -^-J m 



. r 90037575 . 49095 41071 8 

+ [- 2 .6 +-% O^Jm 8 

r_ 1462100355 5463207 9 57165 , 11390625-) , 

~r |_ 2 18 " ' W "2^ + 2 18 — J m 



165044625741 722443913 _ 8198151 705459375 _ 1311187 5 

2 23 + 2 17 2 13 . 5 + 2 20 2 17 



_ 924075-1 
2 15 J 

, r_ 287970294069 87619247043 _ 13794117581 92222296875 
-r|_ 2 212 + 2 20 .5 2 17 .3 2 .5 2 H W 

_ 702232875 _ 34533765 65638125 6476295"! u 

2 19 2 17 2 19 2 lT J m ' 
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By summing the fractions the equation may be written 
mi T 1 . 11 a 2 50625 6 . 822825 r . 65426631 . 

, 514143669 9 . 579596224169 10 , 182494574380633 „ 
+ W ra + ^T5 m + 2*. 3'. 5* m ' 

It will be more suitable for solution if both members are multiplied by 

-*[-!] [l]=4m - m2 -f m5 -« m4 _ 120517 m , _ 4|7J|9 m e. 

The right member then becomes 

50625 7 . 1595025 8 , 112880037 , , 1422559539 10 
—%r~ m + — pr- m 8 H w > m 9 + ^ m 10 

137176160137 „ 47733147493393 12 
+ 2 20 .5 m + 2 22 .3 2 .5 2 m ' 

Calling this K, we have 

c 2 = 1 - i 6* + ,/(0 o -1+| B^f - #, 2 + ^ + i (c 2 - O ) 3 . 
From the preceding expressions for c and 8 , 

i/o" - 0) 3 - - 225 m 3 - 3645 m* - 159429 m' - 1646333 m « • 
* l c ^o; — "23- m — 2T~ m — 2 11 — 2 13 — ' 

whence 
i/i an— H390625 , 553584375 10 60085546875 u 
S (, c — °W — 2>s m 2" m 2 21 m 

1228257320625 ,, 
p 

The substitutions being made in the foregoing value of c 2 and the square root 
extracted, we get 

, , 3 2 201 3 2367 4 111749 , 4095991 , 

c = 1 + m - J2 m 2 - -^ m* - -^ m* - -^ m* - -^-^ m° 

332532037 7 15106211789 8 5975332916861 9 
2i6 32 *" 2 18 3 3 2^ 3 4 

1547804933375567 10 818293211836767367 „ 

— m u — m 

2 25 .3 5 .5 2 28 .3 6 .5 2 ' 

By means of the equation 

1 dio .. 



n dt 1 + m 
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we derive from the last result the ratio of the motion of the perigee to the 
mean motion in longitude, viz. : 

ldw 3 8 . 177 3 . 1659 4 . 85205 5 , 3073531 6 
ndt=T m 2 + ^ m '+ 27- m4 + ^i^ m5 + ^i rF -m« 

, 258767293 . . 12001004273 „ . 4823236506653 „ 
+ 2 18 . 3 2 + 2 IS T3 S m + W7W 

1258410742976387 10 667283922679600927 u 
+ 2 25 . 3* . 5 m + 2 28 . 3 6 . 5 2 m ' 

In this we make the substitution m = = , in which m is the param- 

1 — m r 

eter usually employed, and prolong the resulting series only to the 9th power 

of m. We obtain 

ldw 3 2 . 225 3 . 4071 4 . 265493 s . 12822631 6 
n-3t = * m +-W m + -&- m +-ZF- m + 2". 3 m 

. 1273925965 T 66702631253 8 , 29726828924189 , 
i 2 16 3 2 m 2 18 . 3 s m 2 s 3 4 m ' 

The two terms ending this series are identical with M. Andoyer's and there 
can be no doubt as to their correctness. I do not push this series to the terms 
involving m 10 and m u , as I think the former in terms of the parameter m is to 
be preferred. 

The series which has just been obtained is unsatisfactory on account of 
its slow convergence. It would be of great utility to transform it in such a 
manner that the convergence should be sensibly augmented. Here it seems 
no course is open but to experiment. Confining our attention to parameters 

of the form = , we may seek the value of a which brines about the 

1 — «m ° 

greatest improvement in convergence. It is plain that the adoption of a small 

value for this quantity would not sensibly change the series in this respect, but 

as a is augmented we shall reach a value where one of the numerical coefficients 

vanishes ; if the latter belong to a high power of m, the adjacent coefficients 

will be small. This is true on the assumption that the series tends to become 

a geometrical progression. In the present case it appears that the coefficient 

of m 4 is the first to vanish with augmenting a. Desiring therefore that all the 

coefficients may still be positive after the transformation, I adopt a value of a 

which is less than the value which makes the mentioned coefficient vanish. 

The new parameter adopted is 

m m 



m 



f m 1 — Jm' 
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By making the denominator of a, 4, we secure the advantage that the denom- 
inators of the coefficients of the series are not augmented. With this param- 
eter then, we have the following series : 

1 dco 3 m2 , 141 . , 57 . 41213 . 243353 , ,84226279 m7 

. 1317113479 . , 1125417061277 m9 . 115179069708721 ., 
+ WTW m + 2^73* m + WKW75 m 

, 106545423308527477 mil 

+ 2 28 . 3 9 . 5 2 ' 

The coefficients here diminish more rapidly than in the series proceeding 
according to powers of m. 

Correspondent to the values of n and n' employed in my previous memoirs 
on the lunar theory we have m = 0.0860678013. Substituting this in the 
series just given, we obtain the following result, exhibited term by term ; it has 
been assumed that the series to start from the term involving m 11 may be re- 
garded as a geometrical progression having the ratio £, 

* ^ = 0.0055557498 + 0.0028092554 + 0.0000977435 + 0.0000950403 
n at 

+ 0.0000080501 + 0.0000049959 + 0.0000011207 + 0.0000004292 

[remainder] 

+ 0.0000001260 + 0.0000000418 + 0.0000000209 

= 0.0085725736 . 

The value deduced, without resorting to any expansion in powers of a param- 
eter, is 0.0085725730. The difference of 6 units may be attributed to the 
uncertainty in the estimation of the remainder or to accumulated error in form- 
ing the sum of the terms of the series. 

Oct. 18, 1894. 



